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•/^ ■ Abstract 

We revisit a problem considered by Chow and Hale on the existence of subharmonic solutions for 

rr^ ' perturbed systems. In the analytic setting, under more general (weaker) conditions, we prove their 

results on the existence of bifurcation curves from the nonexistence to the existence of subharmonic 

solutions. In particular our results apply also when one has degeneracy to first order — i.e. when 

^ ' the subharmonic Melnikov function vanishes identically. Moreover we can deal as well with the case 

"ti I in which degeneracy persists to arbitrarily high orders, in the sense that suitable generalisations to 

higher orders of the subharmonic Melnikov function are also identically zero. In general the bifurcation 

curves are not analytic, and even when they are smooth they can form cusps at the origin: we say in 

this case that the curves are degenerate as the corresponding tangent lines coincide. The technique 

we use is completely different from that of Chow and Hale, and it is essentially based on rigorous 

^ ' perturbation theory. 
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O ; 1 Introduction 

O . ,™rT^ 

\^ , Subharmonic bifurcations have been extensively studied in the hterature [TTl |23j. The problem can 

f^ ' be formulated as follows. Consider a two-dimensional autonomous system, and suppose that it has a 

periodic orbit of period T — 2'Kp/q, where p, q are relatively prime integers. Then one can be interested 
in studying whether, under the action of a periodic perturbation of period 27r, some periodic solutions 
exist. Solutions with this property are called subharmonic solutions of order q/p. 

Assume also that the perturbation depends on two parameters. A typical situation is when dissipation 
is present in the system (23 Ell; in this case two parameters naturally arise: the magnitude of the 
perturbation and the damping coefficient. An interesting problem is then to study the region in the space 
of parameters where subharmonic solutions can occur and to determine the bifurcation curves, which 
C^ \ divide the regions of existence and non-existence of these solutions. Such a problem has been considered 

for instance by Chow and Hale |11| . They found that, under suitable assumptions on the unperturbed 
system (essentially a local anisochronicity condition) and on the perturbation, the bifurcation curves exist, 
are smooth and intersect with distinct tangent lines at the origin. The condition on the perturbation, 
if one takes the magnitude of the perturbation as one of the parameters, can be formulated in terms of 
the so-called subharmonic Melnikov function ^29, 24^ . It requires in particular that this function depends 
explicitly on the initial phase Iq of the unperturbed periodic solutions which persist under perturbation. 

In this paper we recover the same result by Chow and Hale, in the analytic setting, and we show that 
the condition on the perturbation can be weakened. In particular the subharmonic Melnikov function can 
be independent of t^. As a consequence the bifurcation curves can be degenerate, in the sense that they 
can have the same tangent at the origin, so that they form a cusp at the origin. Moreover, in general, 
they are not smooth. However if some further assumption is made they turn out to be analytic. 
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In the case of dissipative systems in the presence of forcing, such as those studied by Hale and Taboas 
|25[lll| . our result is significantly stronger as it requires no assumption at all on the periodic perturbation. 
In particular we find the following result. Given any one-dimensional anisochronous mechanical system 
perturbed by a periodic forcing of magnitude e and in the presence of dissipation, there can be analytic 
subharmonic solutions of order q/p only if the dissipation coefhcient 7 is below a threshold value 70 {q/p, e) . 
Here we show that for any rational value p/q there is an integer exponent m — m{q/p) G M* such that 
70(9/^, e) = 0{e'^). This can be related, in a more general context, to a conjecture proposed in [ij. 
Moreover the case m[p/q) = cx) corresponds to infinitely many cancellations, one at each perturbation 
order, which makes such a case very unlikely. Therefore, up to these exceptional cases, we can say that any 
resonant torus with frequency commensurate with the frequency of the forcing term admits subharmonic 
solutions of the corresponding order. In other words, existence of any subharmonic solutions holds without 
making any assumption on the periodic perturbation, other than smoothness. 

Our method is completely different from that of Chow and Hale. It is based on perturbation theory 
More precisely we study the perturbation series of the subharmonic solutions: first we find conditions 
sufficient for these series to be well-defined to all orders, then we prove that if the perturbation is small 
enough convergence of the series can be proved. Technically, this is achieved by using the tree formalism, 
which has been originally introduced by Gallavotti |14j . inspired by a pioneering paper by Eliasson ;l3) . 
and thereafter has been applied in a long series of papers on KAM theory jHl El El El 1211 El HOI 1^ ; 
see also ^H] for a review. We note that with respect to these papers in our case the analysis is much 
easier as we deal with periodic solutions instead of quasi-periodic solutions. In this respect our analysis 
could be considered as a propaedeutic introduction to the tree formalism, in a case in which there is no 
small divisors problem, so that no multiscale analysis has to be introduced; see also JTIJSI for a similar 
situation. In particular Chow and Hale's assumptions on the perturbation reflect a case in which a first 
order analysis is enough to deduce existence of subharmonic solutions. By contrast our results allow the 
analysis of cases in which it can be necessary to go beyond the first order, in principle to arbitrarily high 
orders. 

We also argue that in physical applications it can be essential to have a stronger result. Indeed, in 
a concrete example in which, for instance, the perturbation is a trigonometric polynomial. Chow and 
Hale's assumptions on the perturbation, even if they are generic, fail to be satisfied for most values of 
the periods T. For those values a first order condition is not sufficient to detect the existence of the 
subharmonic solution, and one must go to higher orders. The numerical simulations performed in PP 
for a driven quartic oscillator in the presence of dissipation show that this is necessary if one wants to 
explain the numerical findings for some values of the parameters. 

Our results should be compared also with j^lEl, where a different scenario, such as the persistence 
of the whole invariant manifold corresponding to the resonant torus, arises in a case in which the subhar- 
monic Melnikov function vanishes identically. Our analysis shows that a situation of this kind is highly 
non-generic. 

2 Main results 

Consider the ordinary differential equation 

{a = uj{A)+eF{a,A,C,t), 

\A^eG{a,A,C,t), ^ ' ' 

where (a, ^) e A^ := T x VF, with VK C E. an open set, the map A — + ^{A) is analytic in A, and the 
functions F and G depend analytically on their arguments and are 27r-periodic in a and t. Finally, e, C 
are two real parameters. 

One could also introduce a further (analytic) dependence on e in the functions F and G, and the 



forthcoming analysis could be easily performed with some trivial adaptations. Therefore all the results and 
theorems stated below and in the next sections hold unchanged in that case too. Then, the formulation 
given in |11| is recovered, as a particular case, by introducing the parameter 7 — eC, and setting 
[1. = {fii,^2), with fii = £ and ^2 = 7- 

For e = the variable A is kept fixed at some value Aq, while a rotates with constant angular velocity 
uj^Ao). Hence the motion of the variables {a,A,t) is quasi-periodic, and reduces to a periodic motion 
whenever lu{Aq) becomes commensurate with 1. Define ao(t) = ^(Ao)^ and Ao{t) = Aq: in the extended 
phase-space M x M. the solution {ao{t),AQ{t),t + to) describes an invariant torus, which is uniquely 
determined by the "energy" Aq. If uj{Ao) is rational we say that the torus is resonant. The parameter to 
will be called the initial phase: it fixes the initial datum on the torus. 

As a particular case we can consider that (A, a) are canonical coordinates (action-angle coordinates), 
but the formulation we are giving here is more general. In general all non-resonant tori are completely 
destroyed under perturbation, if no further hypotheses are made on the perturbations F, G (such as that 
the full system is Hamiltonian) . Also the resonant tori disappear, but some remnants are left: indeed 
usually a finite number of periodic orbits, called subharmonic solutions^ lying on the unperturbed torus, 
can survive under perturbation. 

Denote by Tq{A) = 2'it/uj{A) the period of the trajectories on an unperturbed torus, and define 
uj'{A) := dijj{A)/dA. If tj(Ao) = p/q e Q, call T = T{Ao) = 2nq the period of the trajectories in the 
extended phase space. We shall call q/p the order of the corresponding subharmonic solutions. Define 



M{to,C) = ^f dtGiaoit),Ao,C,t + to), (2.2) 



which is called the subharmonic Melnikov function. Here and in the following we do not write explicitly 
the dependence of the subharmonic Melnikov function on ^07 which is fixed once and for all. Note that 
M(io,C') is 27r-periodic in t^. 

We make the following assumptions on the resonant torus with energy Aq. 

Hypothesis 1 One has ijj'{Aq) ^ 0. 

Hypothesis 2 There exists an analytic curve t -^ Co{t) from [0, 27r) to E, such that M(io, C'o(to)) = 
and dM{to, Co{to))/dC ^ for all to € [0, 2tt). 

The function Co(io) is also 27r-periodic in ig- We prove the following result. We prefer to state the 
result in terms of the parameter 7 = eC — instead of C — to make more transparent the relation with 

m- 

Theorem 1 Consider the system Y2.1\) and assume that Hypotheses^ and\^ hold for the resonant torus 
with energy Aq .such that uj{Aq) — p/q. There exist Eq > and two continuous functions 71(e) and 72(2)7 
with 71 (0) = 72(0), 71(e) > 72(e) for e > and 71(e) < 72(e) for e < 0, such that h2.1}) has at least one 
subharmonic solution of order q/p for 72(e) < eC < 71(e) when e G (0, eo) and for 71(e) < eC < 72(e) 
when £ e (— eo,0). 

The situation is depicted in Figure 1, in a case in which the two functions 71 and 72 are smooth. 
The graphs described by the two functions are called the bifurcation curves of the subharmonic solutions: 
they divide the plane into two disjoint sets such that only in one of them there are analytic subharmonic 
solutions. 

In general the functions 71 and 72 are not smooth. However, if some further assumptions are made 
on the subharmonic Melnikov function, smoothness (in fact analyticity) in e can be obtained. Denote by 
Co(<o) and Cg (to) the first and second derivatives of the function Co (to) with respect to to. 







Figure 1: Set of existence (grey region) of subharmonic solutions in the plane (e,7), in a case in which the two 
bifurcation curves e — > 71(e) and e — > 72(e) are smooth and have different tangent lines at the origin. 



Hypothesis 3 IJ tm and Im are the values in [0, 27r) for which the Junction Co(to) attains its minimum 
and its maximum, respectively, then Cq {tm)CQ {t m) 7^ 0. 

The following result holds. 

Theorem 2 Consider the system \2.1]) and assume that Hypotheses Q and hold for the resonant 
torus with energy Aq such that uj{Aq) = p/q. There exist eo > and two functions 71(e) and 72(E), 
analytic for \e\ < Eq, with 7i(0) = 72(0), 71(e) > 72(e) for e > and 71(e) < 72(e) for e < 0, and with 
different tangent lines at the origin, such that \2.1]) has at least one subharmonic solution of order q/p 
for 72(e) < eC < 71(e) when e G (0,eo) and for 71(e) < eC < 72(e) when e £ (— eo,0). 

Theorem is analogous to Theorem 2.1 of ^J, Section 11 — in the analytic setting instead of the 
differentiable one — while Theorem ^ requires fewer hypotheses. In particular it applies when Chow and 
Hale's /ife(a) function vanishes identically. In that case the graphs of the two functions 71 and 72 form a 
cusp at the origin: we refer to this situation as a case of degenerate bifurcation curves, see Figure 2. We 
shall also see in Section0]that in fact, under weaker assumptions than those made in Hypothesis|21 we can 
find smoothness of the bifurcation curves, in the following sense: under suitable assumptions there exist 
two analytic functions 71(e) and 72(e) such that 71(e) = max{7i(e), 72(e)} and 72(e) = min{7i(e), 72(e)} 
for e > 0, and 71(e) = min{7i(e), 72(e)} and 72(e) = max{7i(e), 72(e)} for e < 0. We refer to Hypothesis 
0]and Theorem O in Section^ for a precise formulation of the results. 

We shall see in Section^] — cf. Theorem^] — that for p = 1 one has at least 2g subharmonic solutions 
of order q as far as min{7i(e), 72(e)} < 7 < max{7i(e), 72(e)} and at least q subharmonic solutions of 
order q when (e,7) belongs to one of the bifurcation curves, that is when either 7 = 71(e) or 7 = 72(e). 
This agrees with Chow and Hale's Theorem 2.1 in ^Jj in the cases in which it applies. 

Possible extensions of Chow and Hale's results could be looked for in another direction, such as that 
of relaxing the hypothesis on the unperturbed system. This problem has been studied, for instance, in 

|2ElEZIEni- 

The rest of the paper is organised as follows. Sections and 0] are devoted to the proof of Theorems 
n and 121 More precisely, in Section El we show the existence of a subharmonic solution in the form of a 
formal power series, while in Section^lwe prove the convergence of the series, and we also state Theorem 
121 which provides the aforementioned extensions of Theorem |3 and Theorem 0] on the minimal number 



72(e) 




71 (e) 






Figure 2: In general the bifurcation curves are not smooth. Even when they are smooth they can have the same 
tangent at the origin: in this case we say that the bifurcation curves are degenerate. The grey region in the figure 
represents a case in which the bifurcation curves are smooth and both of them have tangent lines parallel to the 
e-axis. 

of subharmonic solutions of order q. In Section |31 we discuss, as an application of our results, the case of 
a forced one-dimensional system in the presence of dissipation: this will lead to Theorems |31 and |^ which 
extend the results of Hale and Taboas |211 . 

Some final comments, and a comparison with the standard Melnikov theory |29[l24j . are provided in 
Sectional In particular we shall formulate two theorems: Theorem [7| corresponds to the result usually 
discussed in the literature |24| . while Theorem |S1 provides an extension to degenerate situations. 

3 Existence of formal power series for the subharmonic solutions 

We look for subharmonic solutions of H2.1|l which are analytic in e. First, we shall try to find solutions 
in the form of formal power series in e 



ait) = ait,e) - ^ £'=««(*), A(i) = A{t;e) = ^£'=ylW(t), 



(3.1) 



fc=0 



fe=0 



where a^°\t) = uj{Ao)t and A''°\t) = Aq, with uj{Aq) = p/q, and the functions a^''\t) and A^'^^t), 
periodic with period T = 27rp for all fc g IN, are to be determined. We shall see that this will be possible 
provided the parameter C is chosen as a function of e, again in the form of a formal power series in e 



C; = C(e) = ^e'=C('=\ 



(3.2) 



fe=0 



Moreover both the solution {a{t),A{t)) and the constant C will be found to depend on the initial phase 
tQ-. in particular one has C{e) = C(£,io) such that C{e,to + 2n) = C(£,io) and C^°^ = Co(io): and, as 
we shall see, a sufficient condition for formal solvability to hold is that Hypotheses Q] and [21 are satisfied. 
If we introduce the decompositions (|3.1|) and (|3.2I) into (|2.1|l and wc denote with W{t) the Wronskian 
matrix for the unperturbed linearised system, we obtain (cf. 1 for similar computations) 



v(fc) 



it) 



Wit) (|:;) 



Wit) 



dTW-\T) 









(3.3) 



where {a^^\ A^^^) are corrections to the initial conditions, and 

'~^ ^^1 ^^"2 J^^3 

F(''Ht)^[Fia,A,C,t + to)f^:^J2 E ' , , , F{ao{t),Ao,Co,t + to) 

I'l+r2+r3=m 

^ a^'^iHO ■ • ■ a^^^'\t) A('=-i+i)(i) . . . yl('=^-i+'-2)(t) c(fc.i+.2+i) . . . C^^^)^ (3.4) 

fci+...+fcm=fe 

with an analogous definition holding for G^'^^t). Here and henceforth, given a function of several argu- 
ments we are denoting by dk the derivative with respect to the fc-th argument; hence, given the function 
F(a, A,C,t + to) we have diF = dF/da, ^a-F = dF/dA, and ^3^ = dF/dC. Note that by construction 
both F(''^)(i) and G^'''^{t) depend only on the coeflicients a('=')(i), A('=')(f) and C(*=') with /c' < k. 

The Wronskian matrix appearing in (|3.3|l can be written as 

mt)=(; "'^f'). (3.5) 

By using H3.5|l in (|3.3|l we have 

a«(i) = a«+tcc;'(v4o)AW+ / drF(^-i)(r) + cj'(Ao) / dr f At'G^''-^\t'), 

Jq Jo Jo (3-6) 

A'-'^^t) = A'^''^ + f drG'-'^-^^T). 
Jo 

We obtain a periodic solution of period T if, to any order /c S M, one has 



(^Q{k-i)^ :=1. f dTG('=-i)(T) = (3.7) 

^ Jo 



and 

uj'{Ao)A^''^ + {F^''-^^)+u;'{Ao){g^''-^^)^0, g^''-^^{t)^ [ drG^'^-^^T), (3.8) 

JO 

where, given any T-periodic function H we denote by (H) its mean, as done in H3.7|) . 

The parameters a''^^ are left undetermined, and we can fix them arbitrarily, as we have the initial 
phase io which is still a free parameter. For instance we can set a^*'' — for all fc S IN or else we can 
define a^''' — ak{to) for fc S M, with the constants ak{to) to be fixed in the way which turns out to be 
more convenient for computations: we shall see in the next section a reasonable choice. 

Therefore, if equation H3.7|l is satisfied, we have 



a«(t)=y dr(F(^-i)(r)-(F('=-i)))+c.'(Ao) / dr [g^'^-^^t) - (gC^-^^)) 

^(fe)(i)=^(fe)+g('c-i)(t)_ 

with 

uj'[Ao) 
which is well-defined as uj'{Aq) 7^ by Hypothesis ^ 

So, in order to prove the formal solvability of (|2.1|) we have to check whether it is possible to fix the 
parameter C, as a function of e and to, in such a way that H3.7|l follows for all A: > 1. 



(3.9) 



(3.10) 



For k — 1 the condition (|3.7() reads 

(G(°))=M(io,C) = 0, (3.11) 

and we can choose C = Co(io) so that this holds: this is assured by Hypothesis El 
To higher order fc > 1 we can write 

G('=)(a(i), A(i), C, t + to) - d3Giao{t),Ao, Co, t + to) G^^) + r('=)(a(t), A(t), C, t + to), (3.12) 

where the function r('^^(a(t), A{t), C,t + to) depends on the coefHcients C*^'^ ^ of C with k' < k (and on 
the functions a*^'^ ^(t) and A^'' ^(t) with fc' < fc, of course). In other words, in H3.12|l we have extracted 
exphcitly the only term depending on C^'^' . Moreover one has 

(dsGiaoi-), Ao,Co, ■ + to)) ^^J dtd^Giaoit), A^, Co, t + to) = ^^//(to, Co), (3.13) 

and by Hypothesis El one has Dito) := dM{to, Co(to))/9C ^ 0, so that H3.7|l is satisfied provided C'-'^-' is 
chosen as 

C('=) = -;^(r('^K-),^(-),C,-+to))^C,(to). (3.14) 

Therefore we conclude that if we set Co — Co(to) and, for all fc > 1, we choose a'-'"'^ = Q!fc(to), ^''^^ 
according to (|3.1UII and C*-*^^ = Ck{to) according to (|3.14|l . we obtain that in the expansions (|3.1I) the 
coefficients a^'^^t) and A^''\t) are well-defined periodic functions of period T. Of course this does not 
settle the problem of convergence of the series (13. 1|) and H3.2|l . This will be discussed in the next Section. 

4 Convergence of the series for the subharmonic solutions 

Here we shall prove that the formal power series found in Section converge for e small enough, say for 
\s\ < eo for some eo > 0. Then for fixed e G (— £o, £o) we shall find the range allowed for C by computing 
the supremum and the infimum, for to E [0, 27r) of the function to -^ C(e, to). The bifurcation curves will 
be defined in terms of the function C(e, to) — cf. H3.2|l — as 

7i(£)=£ sup C(£,to), 72(e) =£ inf C(£,to). (4.1) 

toe[0,27r) «oe[0,2ir) 

In general the functions (|4.1|1 are not smooth in e. Wc shall return to this at the end of the section. 

To prove convergence of the series (|3.1|l and (|3.2|l it is more convenient to work in Fourier space. First 
of all let us define ut — 2-k/T = 1/q (note that oj 7^ cl'(^o)) and expand 

F(a, A,C,t + to) = Y. e''"e^"^*+*°^ J^^',^(^, C), (4.2) 

so that we can write 
apa2^^aj^^(a;(Ao)t,Ao,Co(to),t + to) = ^e^'^"* ^ e^"'" {iv^r dl^dl-F,,,^,„{A^,Co{h)), (4.3) 

and an analogous expression can be obtained with the function G replacing F. By the analyticity 
assumption on the functions F and G, we have the bounds 



Qr2Qr3 

— j — rFi,„^aa{Ao, Co(to)) 

— j — j-Gj^o, 0-0(^0, Co (to)) 



<^'9?Q?e-''(l'^«l+l'^°l', (4.4) 



for suitable positive constants P, Qi,Q2, k. 
Then, let us write in (|3.1|) 

.C^) (t) = J2 e"'"*al'=) , AC') {t) = Y, ^""""A^^'^ , (4.5) 



a> 



so that 1^1 9|) becomes 



p(fc-i) ^(fe-i) ^(fc-i) 

^W^,^^ +,'(Ao)^, AW = i2^, (4.6) 



for all 1/ 7^ 0, whereas for v — Q one has 

a^ - Mto) - E ^^^ - -'(^o) E #^' 4'-^ - ^^'^ - E ^^^ = -%TT' (4-7) 

with afe(to) so far arbitrary and A^^^ given by (|3.10() . The Fourier coefficients Fy and G!/ can be 
read from (|3.4|l and the analogous expression for G^^''(t). Hence one has 



pw = [na,Ac,i+to)e) = E E E ^ 



.gif^oto 






F.„(Ao,Co(io)) E <^^---«i'!?^^i'!:Vr'---^tt?+:?^C'('='-i+-+^)...C'('="), (4.8) 

fci+...+A;,„ = fe 



and an analogous definition holds for Gv . 
Furthermore one has 

where 

oc 

r(^-)^[r(a,AC,i + to)]f = E E* E 



(4.9) 



ri+r2+r3=m iyop+CTo9+i^i+...+i^,.i+r2 =0 



ri! 



gicroto 



a? 9? 



|^G.„(Ao, Co(io)) E <'^ ■ • ■ <^^ ^tVt^^ ■ • ■ At\t:^^ G('=-+--^) . . . C-C^"), (4.10) 



^ "^ fcl+...+fc,„=fc 

where * means that the term with ri = r2 = and r^ — \ has to be discarded — cf. H3.12(l . 
Therefore we see from the first equation in H4.7|) that it is convenient to fix 

p(fc-i) ^(fe-i) 

a^{t,)^Y.- +^'(^o)E7^ =^ "'^=0, (4.11) 

so that only the functions A^''>{t) have the zeroth Fourier coefficient. 
In particular for fc = 1 we find 



a( 



>'^ = — E e-"*°P,„,.„(Ao,Go(io)) + ^;^^ E e^'^«*°G.„,.„(Ao,Co(io)), 
Ai'^ = — E e"°*°G,„,<,„(Ao,Go(to)), (4.12) 



I'o.a'oGi 



for V ^ Q, and 

for V = Q^ while by writing 



^i'^=-^7^ E e-o*°F.„,.„(Ao,Co(io)), (4.13) 



1^1+1^2=0 i^op+o'og=i^i 
+ E E e'-°*°a2G,„,,„(Ao,Co(io))42^)=Ci(io), (4.14) 

we can express C^^^ in terms of the quantities in H4.12|l . 

In order to study the convergence of the series it is convenient to express ah quantities in terms of 
trees. The strategy is very simple: one iterates the relations l|4.6|) . (|4.7|) and (|4.9|1 . which express the 
coefficients of order k in terms of the coefficients of lower order, until we are left only with the coefficients 
of first order, for which the explicit expressions H4.12|l . (|4.13l) and H4.14|l are at our disposal. 

Trees are defined in the standard way. We briefly recall the basic notations, by referring to |15j for 
an introductory review and further details, and also to 16, 20 for a discussion in similar contexts. 

A tree 9 is defined as a partially ordered set of points, connected by oriented lines. The lines are 
consistently oriented toward a unique point r called the root. The root admits only one entering line called 
the root line. All points except the root are called nodes. Denote with V{0) and L{d) the set of nodes 
and lines in 9, respectively, and with |-/j(0)| and |V^(0)| the number of lines and nodes of 0, respectively. 

If a line £ connects two points Di, O2 and is oriented from 02 to Oi, we say that 02 < Oi and we shall 
write £^2 = ^- We shall say also that £ exits O2 and enters Oi. It can be convenient to imagine that the 
line (. carries an arrow pointing toward the node Oi: the arrow will be thought of as superimposed on the 
line itself. 

More generally we write D2 -< Di if Di is on the path of lines connecting O2 to the root: hence the 
orientation of the lines is opposite to the partial ordering relation -<. Along the path from 02 to Oi all 
arrows point toward Oi. In particular all arrows point toward the root. 

Each line I carries a pair of labels (ft-^, 5t), with hi S {a. A, C} and 5i S {1, 2} such that 5i — 1 for 
hi ^ a. We call h^ and 5^ the component label and the degree label of the line (.^ respectively. Given a 
node D call roi, rc2, and ros the number of lines entering carrying a component label h — a^ h — A^ and 
h = C, respectively. Hence, the values of r^i,r„2,'ro3 are uniquely determined by the component labels 
of the lines entering 0. 

We associate with each node two mode labels i^d , ctd G ^ and we also set for convenience 60 = Se^ . 
With each line £ we associate a further label vi G Z, called the momentum of the line, such that 

i^i = i^£, ^ E ('^"' + '^n') , (4.15) 

mev{e) 

with the constraints that Vi — ii he ~ C and lye j^ ii hi — a. The relation 14.15|l expresses a 
conservation law at each node: the momentum of the line exiting is the sum of the momenta of the 
lines entering plus the mode labels of the node itself. Note that the momentum "flows" through each 
line in the sense of the arrow superimposed on the line. 

The trees with all the labels listed above are called labelled trees. Then given a labelled tree 9 we 



associate with each hne (. a propagator 



1 



1 



D{t^y 



hi = a, A, ui ^ 0, 
hi = A, vi = 0, 
hi = C, vi = 0, 



(4.16) 



and with each node a node factor 



?'tJi!?'D2!?'D3! 



;^'^°*«F,„(Ao,Co(to)), 



a, d„ = 1, 



^ "' ' ' e-°*°G.„(Ao,Cofa)), /^«=a, ^a = 2, 



?'di!?'d2!?'o3! 



iVD= < 



[ll^Oj O2 '-'a „icr„ to 



(4.17) 



rvi^-rx,2^-ro3^- 



rai^.r„2lro3l 



G.„(Ao,Co(io)), h„ 



<5b = 1, 



3'"°*°G,„(Ao,Co(io)), /J«=C, 5„ = 1, 



with the constraint that when h„ — C (and 5o = 1) one has either r„3 > 2 or rd + r„2 > 1- This 
constraint reflects the condition * in H4.10|l . 

Finally we define the value of a tree 6 the number 



Val(6') 



n 5.)( n ^0' 

eeL{0) oev{0) 



(4.18) 



which is a well-defined quantity: indeed all propagators and node factors are bounded quantities. 
Call the order of the tree 9 the number 

k{d) = {£ e 1(9) : hi y^ C} , 



(4.19) 



the total momentum of 6 the momentum i'{9) of the root line, and the total component label of 9 the 
component label h{9) associated to the root line. The number of nodes (and lines) of any tree 9 is related 
to its order k{6) as follows. 

Lemma 1 For any tree 6 one has \L{e)\ = \V{9)\ < 2k{9). 

Proof. The equality |i(^)| = 1^(^)1 is obvious by construction. We prove by induction on k the bounds 

'3fc(0)-2, /i(6')=a,A, 



\V{9)\ < 



3k{9) - 1, h{9) = C. 



(4.20) 



For fc = 1 the bound H4.20|l is trivially satisfied, as a direct check shows: simply compare (|4.12|l to H4.14|l 
with the definition of trees in that case. Assume that the bound holds for all k' < k, and let us show 
that then it holds also for k. Call £q the root line of 9 and Do the node which the root line exits. Call ri. 
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r2, and r^ the number of lines entering Dq with component labels a, A, and C, respectively, and denote 
with 01, ... , 0ri+r2+r3 the subtrees which have those lines as root lines. Then 



ri+r2+r3 



\vm = i+ E i^(^^-^ 



(4.21) 



Then if Iq has component label hi^ G {a, A} we have 

|y(6')| < 1 + 3 (fc - 1) - ra - 2 (ri + ra) < 3fc - 3 < 3fc - 2, (4.22) 

. . . + k{0r^+r2+r3) ~ k — 1, whcrcas if Iq has 



by the inductive hypothesis and by the fact that k{9i 
component label h^g = C we have 



\V{e)\ < 1 + 3fc - ra - 2 (n + ra) < 3fc - 1, 



(4.23) 



by the inductive hypothesis, by the fact that k{9i) + . . . + fc(^-ri+r2+r.-i) = k, and by the constraint that 
either ra > 2 or ri + r2 > 1 in such a case — cf. the comment after (|4.17ll . Therefore the assertion is 
proved. ■ 

Define Qk,v,h the set of all trees of order k{ff) = k, total momentum i/{d) = v, and total component 
label h{9) = h. By collecting together all the definitions given above, one obtains that the Fourier 
coefficients al and Al and the constants Ck can be written for all A; > f in terms of trees as 



a 



(fc) 



(fe) 



9ee 



Al'^) 



J2 Val(0), ly^O, ar=0 

E Val(0), CW= E Val(e). 



(4.24) 



9eek,,..A 



9eek 



The proof of (|4.24l) can be performed by induction; cf. Ts" for details. 

The number of unlabelled trees of order k is bounded by the number of random walks of 2k steps, 
hence by 2^*^ [2S| ■ The sum over all labels except the mode labels and the momenta is bounded again by 
a constant to the power k — simply because all such labels can assume only a finite number of values. 
Finally the sum over the mode labels — which uniquely determine the momenta through the relation 
(|4.f5|l — can be performed by using for each node half the exponential decay factor e"'''^l''°l+'°'°'^ provided 
by the bounds (|4.4|l . The conclusion is that we obtain eventually the bounds 



a, 



i^^ < BiBl^e- 



kl/2 



Ai^^ < BiB^e 



k\v\/2 



C^^^ < BiBt 



(4.25) 



for suitable constants Bi and i?2. This proves the convergence of the series H3.1|l and (|3.2|) for |e| < en, 
with Co small enough. Note that with respect to ^S] here the analysis is much easier as there is no small 
divisors problem. 

The construction described above provides also a useful algorithm which can be implemented numer- 
ically in order to compute the solution to any prescribed accuracy (provided e is small enough). 

Now, we come back to the problem of determining the boundary of the set in the plane (£,7), with 
7 = eC, in which there are subharmonic solutions of order q/p. 

We have to find the solutions of (|4.1|l . that is, solve the equation 



0=^C{e,to) 
dto 



^ C'oito) + eC[ito) + e^C'^ito) 



(4.26) 



where q(to)= dCfc(io)/dio. 
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The function to -^ C{e,to) is analytic in to for all |e| < Eq (for which it is defined and analytic in e), 
so that for fixed e the equation H4.26|l can always be solved. It has at least the two solutions to = Tiis) 
and to = T2(e) corresponding to the absolute minimum and to the absolute maximum, respectively, of 
the function C(e,to). In general these solutions are not smooth in e. This proves Theorem^ 

Suppose now that at the value to such that Co(to) — one has furthermore Cg (to) ^ 0. In that case, 
if To = Toi^) is a solution of l|4.26|l — tq is a point of minimum or maximum for C(e, to) — then tq must 
be analytically close to tp. Hence e -^ To{e) is an analytic function of e, so that also e — ^ Ci{e) and 
£ —* C2{s) are smooth (in fact analytic) in e. Therefore also Theorem |21 follows. 

The last observation suggests how to extend Theorem |5| to obtain smooth bifurcation curves when 
Hypothesis 121 fails to be satisfied. 

Hypothesis 4 There exists fc > 1 such that the functions Cp(to) are identically constant in to for all 
p = 0, . . . , /c — 1. Iftrn and tM are the values in [0, 27r) for which the function Cfc(to) attains its minimum 
and its maximum, respectively, then C'^{tm)C'^(tM) 7^ 0. 

The following result extends Theorem[21 as it deals with the case in which the subharmonic Melnikov 
function does not depend explicitly on to, that is CQ(to) = 0. 

Theorem 3 Consider the system \2.1]) and assume that Hypotheses and 0] hold for the resonant 
torus with energy Aq such that lo{Aq) = p/q- There exist Eq > and two functions 71(e) and 72(2)) 
analytic for \e\ < So, with 71 (0) ~ 72(0) and 71(e) ^ 72(0) for all e 7^ 0, such that the two functions 

.^N fmax{7i(e), 72(e)}, e > 0, fmin{7i(e), 72(e)}, e > 0, 

|min{7i(e),72(e)}, e < 0, ^ [max{7i(e), 72(e)}, e < 0, 

have the same tangent lines at the origin, and Y2. 1]) has at least one subharmonic solution of order q/p 
for 72(e) < eC < 71(e) when e £ (0, eo) and for 71(e) < eC < 72(e) when e G (— eo,0). 

The proof follows the same lines as that of Thcorcm[21 The only difference is that up to order fc — 1 the 
initial phase is left undetermined. In fact to first order one has M{to, C) = M{C) = which fixes C = Cq 
(by Hypothesis |2l, while to orders k' — 2, . . . ,k — l the constants Ck are fixed and are independent of to by 
HypothesisEI Then we can write C(e,to) = £i(e) + £2(e, to), with £i(e) = Co + eCi + . . .+e''"^Cfc_i and 
2^2(e,to) = e'°(Cfc(to) + 0(e)), and from order fc on the constants C^ are fixed as functions of to. Moreover 
equation l|4.26|l reduces to = C^(to) + eC^_|_;^(to) + . . .. Therefore we can reason as in the previous case 
(fc — 0) and we find that Cfe(to) has at least two stationary points to = ti and to = t2, corresponding 
to the minimum point and to the maximum point, respectively. By Hypothesis 0] also £2(e,to) has two 
stationary points at Ti(e) = ti + 0(e) and T2(e) = t2 + 0(e), with Ti(e) and r2(e) analytic in e for e small 
enough. Then we can define 71(e) = C(e, Ti(e)) and 72(e) = C(e, T2(e)): by construction, both 71(e) and 
7i(e) are analytic in e for e small enough. If wc define 71(e) and 72(e) according to (|4.27|) then the proof 
of the theorem is achieved. 

Note that in this case the definition H4.27|l coincides with the general definition 14. 1() for the bifurcation 
curves. Furthermore, if we assume Hypothesis instead of Hypothesis 0] then one has 71(e) = 71(e) and 
71(e) = 71(e), so that also 71(e) and 72(e) are analytic, as stated in Theorem [2 

Finally we note that if the functions Cfe(to) are identically constant in to for all fc G Z+ then one has 
C(e,to) = C(e). In this case the two curves 71(e) and 72(e) coincide, and all values of to are allowed. 
This means that the whole manifold corresponding to the resonant torus persists. On the other hand the 
parameter C must be fixed in a very precise way, as a function of e, and any small deviation from that 
value destroys the torus. This result can be compared with pilll)|. where a similar situation is discussed. 

For (e, 7) inside the set of existence of subharmonic solutions one can investigate how many of them 
exist. For p = 1 the initial phase to varies in the interval [0, 27r(7], where To = 2nq is the period of the 
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unperturbed periodic solution. The function C{e,to) has period 2tt in to; so that it is repeated q times 
in the interval [0, inq]. Hence for any fixed value |£| < eq and any C strictly between the maximum and 
the mininum value attained by the function to —>■ C(e,io) there are at least 2g values ti, i = 2, . . . ,2q, 
such that C = C(e, ti). If C coincides with either its maximum or its minimum then there are at least q 
values ti, i — 2, . . . ,q, such that C = C(e, ti). Therefore we can conclude that, for p = 1, inside the set 
of existence of subharmonic solutions there are at least 2q such solutions, as found in 11 , while on the 
boundary of that set there are q of them. 

We can summarise the discussion above in the following statement. 

Theorem 4 Under the same assumptions of Theorem^ assume p = 1. Take |e| < eo, and for such values 
of e let e —> 71(e) and e —>■ 72(e) be the two bifurcation curves whose existence is assured by Theorem 
Q For min{7i(e), 72(e)} < 7 < max{7i(e), 72(e)} there at least 2q subharmonic solutions of order q. If 
either 7 ~ 71(e) or ^ — 72(e) one has at least q subharmonic solutions of order q. 

Theorem 01 should be compared with Theorem 2.1 in jll| . 

5 Application to dissipative systems with forcing 

Let us consider a one-dimensional system, subject to a conservative force g{x), in the presence of dissi- 
pation and of a periodic forcing. If the periodic forcing and the dissipation coefficient are both small we 
can write the equations for the system as 

x + g{x)'V-ix = £f{x,t), l = eC, (5.1) 

where £f{x, t) is the forcing of period 27r and C is a parameter. Assume that both g and / are analytic 
in their arguments. If / depends only on t, equation (15.1(1 reduces to the equation studied in |25j . 

Let us assume that the unperturbed system (e ^ 0) is Liouville-integrable and anysochronous. This 
means that, in action-angle variables, the equations (|5.1|) can be written in the form H2.1|) . and, further- 
more, that Hypothesis n is satisfied. 

We define the subharmonic Melnikov function in terms of the action-angle variable as in ((2.2() . To 
check that Hypothesis El is also satisfied we use the following result. 

Lemma 2 The subharmonic Melnikov function is invariant under a transformation of coordinates. 

Proof. Consider a system of differential equations in R^ 

x^f{x)+eg{x,t), (5.2) 

and define the subharmonic Melnikov function [23 1241 ITT] for a subharmonic solution xq {t) of period T 
as 

M{to) = ^J dt{fiixo{t))g2ixo{t))~f2{xoit))giixo{t))). (5.3) 

Take the transformation of coordinates ^ ^ x — h{^) . In the new coordinates the system reads 

^^m + niU), (5.4) 

where 0(^) — dh~^{h{^)) f{h{^)) and 7(^) = dh^^ {h{(_)) g{h{^)) , and the subharmonic Melnikov function 
becomes 

M{to)^^l dt(</.i(Co(0) 72(^0(0) -<^2(Co(i)) 71 (^o(t))), (5.5) 
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where ^o(i) is the subhamionic solution expressed in the new variables. 
By noting that 

..-(MO) ^ WO)- 4 (*;■*>) -^Ik?)' '"' 

where J — detdh = 9i/ii92^2 — dih2d2hi is the Jacobian of the transformation, one obtains 

-^(^o) = ^ / dii ( (52/12/1 - 92/11/2) (-5i/i25i + dihm) - 

(-9i/i2/i + dihif2) (92/1231 - 92/1152) 

1 /"^ 1 
= 7^/ dt- {dihid2h2 - dih2d2hi) {fig2 - f2gi) : (5.7) 

where the function /i is computed in £,o(t) and the functions /, g are computed in xo(i) = /i(^o(/))- Hence 
(|5.3|) yields 7VJ(io) = M{to), so that the assertion follows. ■ 

This means that we can compute the subharmonic Melnikov function for the system H5.1|l in the 
coordinates {x,y) = {x,i). In that case the unperturbed vector field is {y, —g{x)) and the perturbation 
reads (0, —eCy + ef{x, t)), so that the subharmonic Melnikov function becomes 

M(io, C)^^J dtyo{t){- Cyoit) + f{xo{t), t + to)) = -C(y2) + (2/o/(xo(-), ' + to))- (5.8) 

Therefore the subharmonic Melnikov function vanishes provided C — Co(/o), where Co(/o) = 
{{yo))~'^ (yof {xo{-) , • + to)), which is well-defined because (j/q) > 0. Moreover one has dM{to,C)/dC = 
~{Vo) 7^ 0. Therefore Hypothesis |21 is also satisfied, and Theorem [3 applies to the system (|5.1f) . 
We can state our result as follows. 

Theorem 5 Consider the system 1^5.1]) and assume that Hypothesis Q holds for the invariant torus with 
energy Ao such that uj{Ao) = p/q. There exist Eq > and two continuous functions 71(e) and 72(e), 
with 71(0) = 72(0), 71(e) > 72(e) for e > and 71(e) < 72(e) for e < 0, such that \2.1]) has at least one 
subharmonic solution of period 2Trp for 72(e) < eC < 71(e) when e £ (0, eo) and for 71(e) < eC < 72(e) 
when e £ (— eo,0). 

Of course Theorem is a corollary of Theorem ^ It should be compared with Corollary 2.3 in JI] 
(cf. also |2ni)- Our result is stronger as it requires, in Chow and Hale's notations, only Hypothesis (Hi), 
which corresponds to our Hypothesis ^ If one assumes also Hypothesis (H4) of JT], which corresponds 
to our hypothesis then Theorem |2 applies, and the result of m] is recovered. 

One expects that, in the case of system (I5.1|l . the two bifurcation curves 71(e) and 72(e) contain the 
real axis, that is min{7i(e), 72(e)} < < max{7i(e), 72(e)}. Indeed for 7 = the equation (|5.1(l describes 
a quasi-integrable Hamiltonian system, and existence of periodic solutions is well known in this case, at 
least under some non-degeneracy condition on the unperturbed system, such as Hypothesis ^ If Co(io) 
is not zero then it is easy to check that the set of existence of subharmonic solutions includes the real 
axis. Indeed, this follows from the following result. 

Lemma 3 The function Co(/o) has zero mean. 

Proof. Call 

F{,xo{t))^j ^f{xo{t),t + to)^J ^f{xo{t),to). (5.9) 
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By (|5.8() the mean (with respect to to) of Co(io) is 

/ lT^Co{to) = -^ -^ -Xo{t)f{xo{t),t + to) 

= £'^ioit)Fix„it)), (5.10) 

which vanishes, as the integrand can be written as a total derivative with respect to t. ■ 

In particular Lemma |31 implies that if Co(io) is not identically constant then its maximum is strictly 
positive and its minimum is strictly negative, hence max{7i(£), 72(e)} > and min{7i(e), 72(e)} < 0. 

To extend the same result to the case in which the functions Ck'{to) are identically constant in to for 
all fc' < /c — 1, with fc > 1 arbitrarily high, is more delicate, and it requires some work. One can reason 
as follows. 

Lemma 4 Assume that for some k G Z+ the coefficients Ck'(to) vanish identically for all k' = 0, . . . ,k — 
1. Then C^(to) has zero mean in to- 



Proof. Write the system (I5.1|) in action- angle variables. Then there exists a Hamiltonian function 
H{a,A,t,£) = Hq{A) +eHi(a,A,t) such that uj{A) = dAHo{A) and 

i a = Lo{A) + edAHi{a, A,C,t) + eC <^{a, A), 

\A = ~edaHi{a,A,C,t)+eC-^{a,A), ^' ' 

where $ = —yda/dy and '^ = ydA/dy. Then (|4.6(l become 

a^'- — idAHt'A +j{Ao)-i-j-dM'~'A +(c$)^^ 

lUJV \ J V (lUJvY \ J V 

A'^^^ ^ — {-d^H^^-^A +{C^f^-^\ (5.12) 

for all fc e M and all ;/ 7^ 0. Moreover (|4.9(l reads 

^ G^.^t'^ + f ('^^ = 0, f ('^^ = (-9„iJj'=-^))^ , (5.13) 

fe'=0 

which, for k = k, gives Fq ' = f q ' and C^^Pg + f q ' = because Ci = . . . — C^_i = by assumption. 
Moreover VP^"^ = — (j/q) ¥" 0, by Lemma|21and Hypothesis El 

Therefore Ck = C'*^^ , with C'-*^-' given by the sum H4.24|l of tree values. We can split the set Ok.o.c into 
the union of disjoint families T as follows. Given a tree G Qk,o,c call Oq the node which is connected 
to the root through the root line, and define Vo{9) as the subset of nodes G V{9) such that all the lines 
£ along the path connecting to Oq have vg ^ 0. Then define T = T{9) as the set of trees obtained from 
9 by "shifting" the root line to any node in Vb(0), i.e. by attaching the root line to any node G Vq{9). 
Of course, as a consequence of the shift of the root line from Oq to D, the arrows of all lines along the 
path between the two nodes are reversed. If one recalls the diagrammatic rules introduced in Section 
0]to associate with any tree 9 a value Val(6'), this means that all lines with labels (/i, (5) = (a, 1) are 
transformed into lines with labels {h, 5) = {A, 1). Moreover the momenta of all such lines change sign. 
The latter property can be seen as follows. The momentum is defined as the sum of all mode labels of 
the nodes preceding the lines — cf. H4.15|l — and the sum of all the mode labels is zero for any tree 
9 G Qkfl,c'- then, when the arrow of a line £ is reversed the nodes preceding i become the nodes following 
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/ ivo+o-o) = =^ ^ (i^d+o'd) = => 2_] f^D = - 2J '''<" (5.15) 



i and vice versa, so that i^£ becomes — i^^ . Hence the propagators of the hnes £ with 61 = 1 change sign, 
whereas the propagators of the Hnes i with 61 = 2 are left unchanged. As a consequence, for each tree 
0' G ^(^) we can write Val(0) = WnVa^fl), where is the node e Vq{9) which the root Une exits and 
Val{9) is the same quantity for aU 6' e ^{d)- Therefore 

Y, Val(6') = Val(6') ^ iiy„. (5.14) 

Moreover one has 

E 

oev(e) DeVo(e) BeVo(e) DGyo(e) 

so that the mean in to of (|5.14|) gives 

"It E Val(0)= r|^VaI(^) 5: ..„ = - r|£val(^) ^ ,,„ = 0, (5.16) 

because the mean is the sum over all labels a„ G ^(^) such that J2vev(e) '^<3 ~ '^vev (e) '^0 ~ ^- ^J 

~ (k) 

using the fact that the set Qk.o.c can be written as a disjoint union of the sets J-', we obtain that Fg 
has zero mean in io, so that the assertion follows. ■ 



Lemma 5 Assume that for some k 6 Z+ the coefficients Ck'{to) are identically constant for all k' — 
0, ...,A:-1. Then Cy {to) = for all k' ^ 0, . . . ,k - 1. 

Proof. The proof is by induction. Fix < k < k, and assume that Ck'{to) = for all k' < k — 1. Then 
by Lemma 0] the function Ck{to) has zero mean. Since it is constant by hypothesis then Ck(to) =0. ■ 

Let k £ Z+ be such that Cfc'(io) is identically constant in to for k' = 0, . . . ,k — 1 whereas Cfe(io) 
depends explicitly on tg- If fc = this simply means that Co(to) depends explicitly on ig. By Lemma |21 
one has Cfe'(tg) = for all fc' < fc — 1, and by Lemma 21 the function Cfc(to) has zero mean in ig- Since 
Ck{to) is not identically constant then supj^^grg 27r) C'fc(^o) > and infj^jgjg 27r) Ck{to) < 0. Furthermore, 
in such a case C(e,io) = £'^{Ck{to) + 0(e)), so that also 

sup C{s,to)>0, inf C{s,to)<0, (5.17) 

toe[0,2Ti-) toe[0,27r) 

for e small enough. If we recall the definition (|4.1|) of the bifurcation curves we can formulate the following 
result. 

Theorem 6 Under the same assumptions of Theorem\^let e — > 71(e) and e —^ 72(e) be the two bifur- 
cation curves whose existence is assured by Theorem\^ One has 71(e) > > 72(e) for e G (0,eo) and 
7i(£)<0<72(e)/oree(0,eo). 

As (|5.17ll shows, if there is fc > such that Cfe'(to) = ioi k' = 0, . . . , k — 1 and Ck(to) 7^ 0, then one 
has the strict inequalities 71(e) > > 72(e) for e G (0,eo) and 71(e) < < 72(e) for e G (0,eo). On the 
contrary if all Ck vanish identically, so that the full function C(e, ta) has to be zero, then 71(e) = 72(e) = 0. 

Therefore Theorems |31 and El show that any one-dimensional anisochronous mechanical system, when 
perturbed by a periodic forcing and in the presence of dissipation, up to the exceptional cases in which the 
functions Ck{to) are constant — hence vanish, by Lemma[Sl — in tg for all k G Z^, admits subharmonic 
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solutions of all orders, without any assumption on the perturbation, — a result which does not follow 
from the analysis of |25l ITT] . 

The case that all the functions Cfe(io) are identically constant in Iq is really exceptional. This can 
be appreciated by the following argument. If the function C(e,io) does not depend on ip then not 
only, by Lemma |S1 it must vanish identically, i.e. C(e,io) = C'(e) = 0, but we find also that Iq is left 
undetermined. In other words the periodic solution persists for all values of io- This means that if we 
take the system (|5.1(l with 7 = 0, so that it becomes an autonomous quasi-integrable Hamiltonian system, 
with no dissipation left, the full resonant torus with frequency lu — p/q persists under perturbation. This 
situation is certainly unlikely — even if not impossible in principle. For instance one can take the system 
described by the Hamiltonian 
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H{x,y,t)^-y' + ^x^ + efit)i^-y' + -x''-E^ , (5.18) 

with E corresponding to the unperturbed solution (xo(t),yQ{t)) with frequency w. Then such a solution 
still satisfies the corresponding Hamilton equations for all values of e and for all values of the initial phase 
to', that is the full resonant torus with frequency uj persists. In particular ii lo — p/q is rational — so that 
the frequency of the unperturbed solution becomes commensurate with the frequency 1 of the perturbing 
potential /, the corresponding torus is resonant. 

It is important to stress that if we look for a subharmonic solution which continues some unper- 
turbed periodic solution with a given period T = 2'Kq/p it can happen that the corresponding integral 
{Uofi^oi'):' + ^0)) identically vanishes. In fact, if / is a trigonometric polynomial (which is often the 
case in physical applications) this happens for all p/q but a finite set of values. An explicit example has 
been considered in fTl. In these cases the subharmonic Melnikov function does not depend on to and it 
is linear in C: hence (|5.8|l can be satisfied only by taking Co (to) = 0. Then, it becomes essential to go 
to higher orders of perturbation theory to study for which values of C a subharmonic solution of order 
q/p appears. Again, we refer to ^ for a situation in which one must perform a higher order analysis to 
explain the numerical findings. 

6 Conclusions, and final comments 

The Melnikov theory [5^ considers systems which, in suitable coordinates, can be written as in p. 1(1 . 
without the parameter C: 

ia = ujiA)+eF{a,A,t), 

\A^eG{a,A,t), ^ ' ' 

where all notations are as explained after (|2.1|) . Define the subharmonic Melnikov function as 

M{to)^l-l dtG{ao{t),Ao(t),t + to), (6.2) 



and set M'(io) = dM{to)/dto. 

We can repeat the analysis of formal solvability in Section O with some adaptations due to the fact 
that no extra parameters C'*'' are at our disposal to any perturbation orders. 

In particular to first order one needs M(io) = 0, so that to must be a zero for the subharmonic 
Melnikov function. To higher orders we can write 

G(^)(a(i), A(t), C, t + to) = diG{aoit),Ao,t + to) a^^^ + T^''\a{t), A{t), t + to), (6.3) 

where the function F*^'') (a(t), A{t),t + to) depends on the corrections a'-'^ •* to the initial phase, only with 
k' <k. 
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To any perturbation order k the constant a^''' is left undetermined. Anyway we are no longer free to 
fix it equal to some arbitrary value, for instance zero, as we no longer have the initial phase to and the 
constants C'^'^-' as free parameters. Hence we shall need the corrections a^'^-' to assure solvability of the 
equations of motion to any order. This will be possible in the light of the following results. 

Lemma 6 One has uj{An){diG{ao{-), Aq, ■ + to)) = -Af'(to). 

Proof. One has 

—G{ao{t),Ao, t + to) = Lo{Ao) diG{ao{t), Ao, t + to) + --G{ao{t), ^o, t + to), (6.4) 

at oto 

where we have used the fact that Ao{t) = and ao{t) — ^(^o)- If we integrate Ht).4|l over a period we 
obtain 

0=^ / dt^G(aoit),Ao,t + to)=io{Ao){diGiaoi-),Aor+to)) + ^{G{aoi-),Ao,- + to)), (6.5) 
T Jo at dto 



so that 



u;iAQ)diGiaoit), AQ,Go,t + to) = - J-(G(ao(i), ^o, t + to)) - -M'(io). (6.6) 

oto 



Hence the assertion follows. ■ 

Thus, if we impose the condition that to be a simple zero for the subharmonic Melnikov function we 
find that in (|6.3|) the derivative diG{ao{t),Ao,t + to) is different from zero, and this allows us to fix a^''^ 
in such a way as to make the mean of G^''\a{t), A(t),t + to) vanish. Hence by fixing the constants A^''^ 
as explained in Section|31and the constants a^*^-* as stated above we find that a solution in the form of a 
formal power series in e exists. The convergence of the series, hence the existence of an analytic solution, 
can be proved by reasoning as in Section 0] We do not repeat the analysis, which would essentially be a 
word for word copy of what was done in Section ^ 

Therefore we have proved the following result — well-known in the literature |24| . 

Theorem 7 Consider a periodic solution with frequency lu — p/q for the system f6.1p . and assume that 
to is a simple zero for the subharmonic Melnikov function 1^6. S^) corresponding to such a solution. There 
exists £o > such that for \s\ < So the system i6.1\} has at least one subharmonic solution of order q/p. 

However, our analysis permits us to generalise the result above. Define 

Afo(to) = M(to), M,(to) = (rW(a(.),A(.),-+io)), fc e M, (6.7) 

where the notations of 16.31) have been used. Note that if Mk'{to) vanishes identically for all k' — 
0, 1, . . . , A; — 1, then Mk{to) is well-defined. The following result follows. 

Theorem 8 Consider a periodic solution with frequency lu ~ p/q for the system i)6'.J|) . Assume that the 
functions Mk' are identically zero for all k' — 0,1, ... ,k — I, and assume that to is a simple zero for the 
function Mk{to). There exists eo > such that for \e\ < eo the system 16'. J |) has at least one subharmonic 
solution of order q/p. 

Of course also the system (|2.1|l can be studied as illustrated in this section. One simply treats the 
parameter C as fixed, and one fixes the initial phase to in such a way that Theorem [7| or Theorem |S1 can 
be applied — of course, provided the corresponding hypotheses are satisfied. This has been done in PP 
to study the subharmonic solutions of a forced cubic oscillator in the presence of dissipation. 
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We also note that, as a particular case of Theorem|Sl it can happen that Mk{to) = for all k G Z+. In 
that case formal solvability of the equations holds to all orders, and the convergence of the series requires 
no condition on Iq, and it can be proved by proceeding as in Sectional In particular in such a case the 
full resonant torus persists under perturbation. Of course, the vanishing identically of all functions Mk 
is a very unlikely situation, and, without any further parameter at our disposal, we can hardly expect 
this ever to happen. This shows that the persistence of the full torus when the subharmonic Mclnikov 
function is identically zero is a very rare event. 

The bifurcation curves studied in this paper concern subharmonic solutions which are analytic in e. 
In principle our results do not exclude existence of other subharmonic solutions which are not analytic. 
Indeed, one could wonder whether other periodic solutions with the same period exist for e ^ 0. In 
the presence of dissipation, it is unlikely that solutions other than the attractive ones found with the 
method we have used, would be relevant for the dynamics — cf. for instance the problems investigated in 
01 13 El n El ■ In general the situation can be delicate; for instance when one investigates quasi-periodic 
solutions corresponding to lower-dimensional tori of quasi-integrable systems, where uniqueness becomes 
a subtle problem - cf. for instance |22[I12| . Despite this, there are cases in which the problem can be 
settled -cf. |3|in|. 
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